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Abstract 

Transverse momentum event -by-event fluctuations are studied within the string- 
hadronic model of high energy nuclear collisions, LUCIAE. Data on non-statistical 
Pt fluctuations in p+p interactions are reproduced. Fluctuations of similar mag- 
nitude are predicted for nucleus-nucleus collisions, in contradiction to the prelimi- 
nary NA49 results. The introduction of a string clustering mechanism (Firecracker 
Model) leads to a further, significant increase of pt fluctuations for nucleus-nucleus 
collisions. Secondary hadronic interactions, as implemented in LUCIAE, cause only 
a small reduction of pt fluctuations. 
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1 Introduction 



Experimental results on event-by-event fluctuations in nuclear collisions at relativistic 
energy may serve as a crucial test of various theoretical approaches to the collision process. 
Until now the existing models were developed with the aim to reproduce the data on single 
particle yields and two particle correlations. To which extend can the same models predict 
event-by-event fluctuations? 

A method to study fluctuations of global kinematical observables in high energy nu- 
clear collisions (like fluctuations of total transverse momentum) was proposed in [[[[]. It 
was applied to analyze experimental data || [3| and test theoretical approaches ||, ||. || . 
The method is based on the fact that in elementary interactions (e.g. p+p interactions) 
particles are produced in a correlated way which leads to the observation of large (non- 
statistical) event-by-event fluctuations. These 'elementary' fluctuations provide a scale 
relative to which the fluctuations in nuclear collisions can be studied in a model indepen- 
dent way. The method also provides a statistical tool to account for 'trivial' geometrical 
fluctuations. 

Recently it was shown that initial state scattering models lead to increase of pt fluc- 
tuations M. This is in contradiction with the experimental results which indicate that 
the properly normalized px fluctuations in central Pb+Pb collisions at 158 A-GeV are 
significantly smaller than the corresponding fluctuations in p+p interactions ]2], |7|. To 
which level can the experimental result be understood as an effect of equilibration due to 
hadronic cascading? What is the role of possible collective effects at the early stage of the 
collision, like clustering of strings? The aim of this paper is to investigate these questions 
using the string-hadronic model LUCIAE ||. 

The paper is organized as follows. The analysis method used further to study event- 
by-fluctuations is recalled in Section 2. In the Section 3 we briefly sketch the basic physics 
ingredients of the LUCIAE. Model calculations concerning p? fluctuations are presented 
and discussed in Section 4. Summary and discussion close the paper. 

2 A Measure of Event-by- Event Fluctuations 

Event -by-event fluctuations in nuclear collisions are usually dominated by the trivial 
variation in impact parameter from event to event and by the purely statistical (here we 
mean statistics for classical particles) variation of the measured quantities. An analysis 
method that allows to remove these trivial contributions and to determine the remaining 
part of event-by-event fluctuations of transverse momentum has been proposed in U. 
Following this reference we define for every particle i in an event: 

Zi = Pt x - Pt, 

where pr is the mean transverse momentum of accepted particles averaged over all events 
(the inclusive mean). Using z% we calculate for every event 

N 

1=1 
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where N is the number of accepted particles in the event. With this definition one obtains 
the following measure of event-by-event fluctuations: 



■±- PT 



{Z2) z\ (1) 



\ (N) 



where (N) and (Z 2 ) are averages over all events and the second term in the r.h.s. is the 
square root of the second moment of the inclusive transverse momentum distribution. 
The physical motivation for studying <£> PT was given in Q: experimental data on N+N 
interactions show that particles in these collisions are not produced independently ||. 
One observes large scale correlations that lead to, e.g., a correlation between the event 
multiplicity and the average of the particles. The correlated particle emission in ele- 
mentary processes can be used to probe the dynamics of nuclear collisions by measuring 
to which degree this correlation is changed when going to p+A and/or A+A collisions. 
For this purpose, $ PT as a measure of fluctuations has two important properties. For 
a large system (i.e. an A+A collision) that is a superposition of many independent ele- 
mentary systems (i.e. N+N interactions), <3? PT has a constant value that is identical to 
that of the elementary system. In other words if the central Pb+Pb collisions were an 
independent superposition of N+N interactions, the value of <3> PT would remain constant, 
independent of the number of superimposed elementary interactions in a single event and 
its distribution in the studied sample of the events. If on the other hand the large system 
consists of particles that have been emitted independently, $ PT assumes a value of zero. 
Thus $ PT provides us with a scale characterising the fluctuations in nuclear collisions 
relative to elementary interactions at the same energy. 

One should expect that $ PT is sensitive to both event-by-event fluctuations in the creation 
of the early state of the collision as well as in its subsequent evolution until freeze-out. 



3 The string— hadronic model — LUCIAE 



The LUCIAE model is developed based on FRITIOF |TJ. In the FRITIOF model a 



hadron is assumed to behave like a massless relativistic string corresponding to a confined 
color force field of a vortex line character embedded in a type II color superconducting 
vacuum. In FRITIOF, during the collision two hadrons are excited due to longitudinal 
momentum transfers and/or Rutherford Parton Scattering (RPS). The highly excited 
states will emit bremsstrahlung gluons according to the soft radiation model. They are 
afterwards treated as excitations, i.e. the Lund Strings, and allowed to decay into final 
state hadrons according to the Lund fragmentation scheme [O . 

The LUCIAE includes all elements of the FRITIOF model and, additionaly, two com- 
ponents: a 'Firecracker' model and a model of hadron rescattering. In the Firecracker 
model it is assumed that groups of neighbouring strings may form interacting quantum 
states (clusters) so that both the emission of gluonic bremsstrahlung as well as the frag- 
mentation properties can be affected by the large common energy density of the interacting 
strings. The maximum transverse momentum of the emitted gluons is found to fulfill a 
condition |12[ k± max < \/ fiMtot, where M to t is the total excitation energy of a cluster and 
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fiMtot effectively corresponds to an energy density over a region of transverse size 
Consequently, when it comes to heavy ion collision predictions the Firecracker model will 
correspond to an essential enhancement of (mini)jets in the center of phase space, which 
contributes to high px enhancement [12| . 

In the rescattering model, the produced particles (which consist of hadrons after strong 
decays and the participant nucleons) are distributed randomly in the geometrical overlap- 
ping region between the projectile and the target nuclei. The target (projectile) spectator 
nucleons are distributed randomly outside the overlapping region and inside the target 
(projectile) sphere. A formation time is given to each particle and a particle starts to scat- 
ter with others after it is "formed" . Two particles will collide if their minimum distance 
dmin < J^tot/^, where a to t is the total cross section in fm 2 and the minimum distance is 



calculated in the cms frame of the two colliding particles (for details see ||). 
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4 The p T Fluctuations in LUCIAE 



The analysis of the event-by-event fluctuations in LUCIAE is done in several steps. In 
the first part of this section we study event-by-event fluctuations in p+p interactions. In 
the following subsections the pt fluctuations in A+A collisions are analyzed. We study 
effects of independent string superposition, string clustering and finally the role of the 
hadronic rescattering. 

In the analysis presented in this paper the <ft PT variable was calculated using all charged 
particles stable with respect to strong interaction, charged spectator fragments were ex- 
cluded. The number of generated events for all analyzed event samples was larger than 
10 4 . The statistical error of the <p PT variable was calculated according to formula given in 
the Appendix. 



4.1 Fluctuations in p+p Interactions 

It is well established experimenatlly || that particle production in p+p interactions at 
high energy is correlated. It was pointed out in Ref. [ l] that this correlation could lead 
to a large non-statistical event-by-event fluctuations which can be quantified by the $ 
measurement. 

It is essential that the string-hadronic models which are used for the analysis of 
the event-by-event fluctuations in A+A collisions are first checked as to whether they 
reproduce the corresponding fluctuations measured in p+p interactions. This is because 
in these models the latter process is used as an input for the calculation of the properties 
of the A+A collisions. 

It was shown that the $p T variable is sensitive to the correlation between the form 
of the pt distribution and the event multiplicity. Therefore in Fig. 1 we compare the 
experimental data on (pt) vs n dependence for p+p interactions at 200 GeV |13| with the 



LUCIAE results. The comparison is done separately for 7r + and ir~ mesons. Two versions 
of LUCIAE are used. In Figs, la and lc the results for the full version of the model, 
which includes hard processes (RPS), are shown. The results obtained by switching off 
hard processes in the model are presented in Figs, lb and Id. We note here that the 
string clustering model as well as hadronic rescattering play no role for p+p interactions. 

The comparison indicates that LUCIAE underpredicts both the absolute magnitude 
and the strength of the (pt) vs n dependence. The hard scattering processes, which 
occur in about 5% of p+p interactions at 200 GeV, play a minor role when the (j>t) vs n 
dependence is considered. 

The fact that the (pr) vs n dependence is underpredicted by LUCIAE may suggest 
that the event-by-event px fluctuations are smaller in the model than in the data. This is 
however not the case. The $ PT value calculated within LUCIAE without hard processes 
is about 15 MeV/c for p+p interaction at 158 GeV. This value of $ PT agrees with the 
corresponding preliminary experimental value obtained by the NA49 Collaboration [7j 
when the acceptance effects are taken into account. The introduction of hard scattering 
increases the $p T value to 30 MeV/c leading to an overestimation of the experimental 
number. 

The above results illustrate two important features of $ PT . Its value is determined not 
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only by the (pt) vs n dependence but also by other correlations present in the particle 
production process (e.g. jet production). It is sensitive to different properties of the 
production process than the inclusive or semi-inclusive observables and therefore yields 
additional information otherwise not available. 

The strong increase of the $ PT value when the hard processes are included may be 
understood as due to the fact that the final state particles originating from the hard 
process (jets) are strongly correlated in momentum space. 

In order to trace the origin of the particle correlation in the p+p interactions for soft 
processes in LUCIAE we exchanged the string fragmentation scheme by an independent 
fragmentation scheme in which energy-momentum and quantum charges are not con- 
served. This results in a reduction of the $ PT value by a factor of about 5, to almost zero. 
Thus we conclude that the conservation laws are responsible for large event-by-event 
fluctuations in soft p+p interactions at SPS energies in LUCIAE. 

For the further study of pt fluctuations in A+A collisions we selected the LUCIAE 
version without hard processes as only this version reproduces correctly the magnitude of 
the pt fluctuations measured for p+p interactions |7J . 



4.2 Superposition of Strings 

The fundamental assumption of the string-hadronic models of A+A collisions is that the 
basic physics of these collisions can be pictured as an (almost) independent superposition 
of nucleon-nucleon (N+N) interactions. 

Thus in order to study properties of this minimal model with respect to pt fluctuations 
we analyze the fluctuations calculated within the LUCIAE without hard processes, string 
clustering and hadronic rescattering. For this version of the model we calculate the $ PT 
for central (6 = fm) A+A collisions and plot it as a function of A in Fig. 2. 

One observes a independence of $ PT of A. This behaviour is expected, by definition, 
for a model in which A+A is assumed to be an independent superposition of N+N interac- 
tions. In the studied version of LUCIAE there are two effects which can cause deviations 
from the independent superposition picture. The nucleons in the nucleus have Fermi mo- 
tion and the string excitation increases with the size of the colliding nuclei. Our numerical 
results show that both effects play a minor role when the $ PT is studied. Therefore, in 
this respect, the minimal version of the string-hadronic model of A+A collisions can be 
considered as an independent superposition of N+N interactions. 



4.3 String Clustering 

In the high density stage of the A+A collisions the picture of independent string superpo- 
sition is obviously unrealistic. A modification of this picture is modeled, in the framework 
of string-hadronic models, by introduction of string clustering (Firecracker model) || or 
string fusion [O]. 

In the Firecracker model clustering of strings allows for a collective gluon emission 
and therefore is expected to increase the pt fluctuations. The probability that the string 
will form a cluster increases with the size of the colliding nuclei and consequently the pt 
fluctuations due to string clustering should increase with the size of the colliding nuclei. In 
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fact this is observed in Fig. 2, where the results of the calculations with string clustering 
effect are shown. For central Pb+Pb collisions the value of <E> PT increases by a factor of 
about 3 when the string clustering is added to the minimal version of the LUCIAE model. 

4.4 Rescattering 

The early stage of A+A collisions is followed by the stage in which hadrons and hadronic 
resonances are effective degrees of freedom. Considerable rescattering between them is 
expected to take place before the final decoupling of produced particles. For the inter- 
pretation of the measured event-by-event fluctuations it is crucial to understand the role 
played by the hadronic rescattering process. 

The $ PT calculated for central S+S and Pb+Pb collisions within the minimal version 
of the LUCIAE supplemented by the hadronic rescattering model is shown in Fig. 2. 
One observes that the pr fluctuations are slightly reduced by the hadronic secondary 
interactions. This trend can be, in fact, understood as the rescattering should lead to the 
equilibration of the final state and therefore it should reduce the fluctuations established 
at the early stage. From that point of view it may even be surprising that the role 
of the rescattering is relatively small, as each produced hadron rescatters in average 6 
timesi until decoupling in central Pb+Pb collision modeled by the present version of the 
LUCIAE. 
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5 Summary and Discussion 



The analysis of event-by-event fluctuations in string-hadronic model LUCIAE presented 
in this paper was triggered by the experimental results of the NA49 Collaboration 0] . 
They indicate that pr fluctuations, as measured by $ PT measure, are significantly reduced 
in central Pb+Pb collisions at 158 A-GeV in comparison to p+p interactions at the same 
energy. 

Our main results can be summarized as follows. 

1. The $ PT value in p+p interactions modeled by LUCIAE is greater than zero. This 
is caused by correlations between particles introduced by conservation laws and by 
hard scattering. 

2. The $ PT value is independent of the size of the colliding nuclei when A+A collisions 
are modeled by the minimal version of the LUCIAE model (without hard scattering, 
string clustering and hadronic rescattering) . It remains constant at the value found 
for p+p interactions. 

3. A string clustering effect, as introduced in LUCIAE, causes a strong increase of $ PT 
with increasing size of the colliding nuclei. 

4. Final state hadronic rescattering, as modeled in LUCIAE, only weakly decreases 
the value of <J> PT . 

Our results should be confronted with the recently presented results on event-by-event 
fluctuations as measured by the $ PT . 

It was found in J| that initial state scattering models predict an increase of $ PT with 
the size of the colliding nuclei. This behaviour is similar to the behaviour obtained by 
us for the string clustering effect. The relative contribution of both effects (initial state 
scattering and string clustering) increases with the size of the colliding nuclei. 

The influence of hadronic rescattering on event-by-event fluctuations was studied in 
Ref. using the VENUS model |[J and in Ref. § using the UrQMD model [0]. The 



secondary interactions as modelled by VENUS were shown not to influence fluctuations 
in pseudorapidity measured by & v |J. UrQMD finds a strong reduction of the value 
of $ PT when going from p+p interactions to central Pb+Pb collisions. This result is 
in contradiction with our finding (see point 4 above). We are aware of the following 
differences between the current analysis and the analysis done within UrQMD model 
||, [18]. In the UrQMD analysis all particles at midrapidity (—0.5 < y* < 0.5) are used 



whereas in our analysis only charged particles but without rapidity selection are included. 
We checked however that our conclusion on the weak influcence of rescattering remains 
unchanged when a 'UrQMD acceptance' is used for LUCIAE events. The number of 
central Pb+Pb events with rescattering was 10000 for LUCIAE and 1600 for UrQMDf]. 
The rescattering prescriptions are different in UrQMD and LUCIAE. We checked however 
that the number of secondary hadronic collisions in central Pb+Pb collision is similar 
(about 6 rescatterings per final state hadron). 

5 In both models generation of a single Pb+Pb event takes about 30' CPU time on Pentium Pro 200. 
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The difference in the conclusions may be, also, due to the fact that in UrQMD model 
Pauli blocking of baryons is taken into account. This leads us to a comment on the results 
obtained recently by Mrowczyhski |J. He calculated the value of $ PT for an equilibrium 
ideal gas in the grand canonical approximation. In the case of classical particles $ PT = 
0. The $ PT value for fermions is large negative, whereas for bosons large positive. He 
estimates that due to the dominance of pions in the final state the value of $ PT in the 
hadronic gas should be large and positive and therefore it is not easy to understand the 
NA49 results even assuming full equilibration of the matter. 

We conclude from the above sketchy review that the basic discrepancy between the 
results of various models remains to be reconciled. Vis a vis the preliminary data ob- 
tained by NA49 0, [F| which indicate a vanishing $ PT measure both the results of our 
microscopic string-hadronic model LUCIAE, and expectation based on an equilibrium 
hadron quantum gas miss the mark. On the other hand the UrQMD model which is 
based on a similar physical picture predicts a vanishing $ PT measure. Further, ongoing 
work is devoted to pinpoint the origin(s) of this apparent discrepancy. We also urgently 
await final affirmation of the preliminary NA49 data. 
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Appendix 



A. The fluctuation measure & x , where x is any single particle variable, defined by 
Eq. 1 can be expressed in an explicit way by single event variables N, X and X 2 , where 
N is the number of particles and 

N N 

x = ^Xi, x 2 = J2 x l 
i=i i=i 

It is easy to show that the definition (Eq. 1) is equivalent to 

= /(#»> _ 2(x)(xn) (N*)(x)* y/* _ ((xi _ (xyy/ 2 

x \ (TV) (N) 2 (iV) 3 J \ (N) (N) 2 ) • 1 ' 

The above formula allows to calculate value of $ x during a single pass of data processing, 
without initial evaluation of the inclusive mean, x. 

It is common to study event-by-event fluctuations by the analysis of the variance 
of the ratio X/N (the mean x for particles from a given event). From the definition of 
variance we get: 

We make a trivial observation that the $ x and V(X/N) are different functions of different 
moments of basic single event observables X and N. Thus analysis of & x is not equivalent 
to the analysis of V(X/N). 

B. One can estimate a statistical error of the $ x , a(Q x ), when the single event variables 
N, Z 2 and Z 2 (Z 2 = J2iLi Z T) are considered as the original random variables. The <f> x 
can be then written as: 



*-= 3£n£ - ra^r . (4) 



where iV^y is the number of events. The a($ x ) can be therefore expressed as: 

^ ^ (^)V.^ + (^)V(*,« + (*)V,») + 



<5Z 2 <5iV v ' y 5Z 2 5Z 2 y " ' 5N 5Z L 

where 5<& X /S(N, Z 2 , Z 2 ) is a derivate of $ s over a single event variable (N, Z 2 , Z 2 ) and 
^(Yi,!^) — ({Yi — (Yi))(Y 2 — (Y 2 ))). Of course finally the expression for the statistical error 
of $ x can be written in terms of the algebraic moments as in the case of the observable 
<$> x itself (Eq. 2). 
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Figure 1: The dependence of the average transverse momentum of 7r~ (Figs, la and lb) 
and 7r + (Figs, lc and Id) mesons on multiplicity of negatively charged hadrons for p+p 
interactions at 200 GeV. The experimental results ]13[ are indicated by the filled squares, 
whereas the results obtained within the LUCIAE model by the open diamonds. The hard 
processes are included in the case of Figs, la and lc, and excluded in the case of Figs, 
lb and Id. 
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Figure 2: The dependence of the $ PT on the nuclear mass number of the colliding nuclei 
calculated for A+A collisions at 158 A GeV at zero impact parameter using LUCIAE 
model. Results obtained without hard scattering, string clustering and secondary hadronic 
interactions (minimal version of the model) are indicated by open diamonds. Behaviour 
obtained by including string clustering is shown by filled squares. Filled circles indicate 
the value of $ PT obtained when the hadronic rescattering is added to the minimal version 
of the model. 
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